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We give a quantum benchmark for teleportation and quantum storage experiments suited for pure and mixed 
test states. Tiie benchmarlc is based on the average fidelity over a family of phase-covariant states and certifies 
that an experiment can not be emulated by a classical setup, i.e., by a measure-and-prepare scheme. We give an 
analytical solution for qubits, which shows important differences with standard state estimation approach, and 
compute the value of the benchmark for coherent and squeezed states, both pure and mixed. 
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Introduction: A central question in quantum information 
theory is whether a particular quantum protocol can be real- 
ized with the same efficiency by classical means; ultimately, 
quantum information stands on the advantage of quantum over 
classical systems in performing certain tasks. This question is 
relevant in experimental implementations of such protocols, 
as they become imperfect in real, unavoidably noisy, experi- 
ments. It is then essential to assess whether the same exper- 
imental result could have been obtained by using only clas- 
sical (less costly) resources. So, for example, in a particular 
teleportation experiment — which involves generating entan- 
glement, performing complicated Bell measurements, fighting 
decoherence, etc. — one may ask whether the same goal, i.e., 
mapping the state of a system onto a second system at a dif- 
ferent space-time location, could have been achieved by mea- 
suring the quantum state of the first system, transmitting the 
collected information, and preparing the state of the second 
accordingly. In this letter, we propose and calculate quantum 
benchmarks that certify that a certain implementation of tele- 
portation cannot be realized classically, or more precisely, by 
a measure-and-prepare strategy. Of course, these benchmarks 
also apply to any protocol that can be understood as a realiza- 
tion of the identity channel, e.g., quantum storage. Our bench- 
marks are based on phase-covariant families of states, and thus 
they are computationally manageable and experimentally fea- 
sible (shifting phases is straightforward in, e.g., optical ex- 
periments). Last but not least, our benchmarks apply when 
test states are mixed, which is the standard situation in exper- 
iments. Previously proposed quantum benchmarks Uj, |2|, |3|] 
were either restricted to pure states or were not strict bounds, 
since they were based on sub-optimal classical strategies (see 
footenote 15 in [30). We give rigorous quantum benchmarks 
for both pure and mixed Gaussian states that can be imme- 
diately applied to current experiments on continuos-variable 
(CV) systems, such as optical fields and atomic ensembles 
ii. 

Background and methods: The literature on quantum 
benchmarks originated in the context of CV teleportation 
experiments, the first of which was performed in 1998 by 
A. Furusawa et al. 14(1. In this experiment, optical coherent 
states were teleported using squeezed-state entanglement. As 
benchmark, they used the average fidelity, T, that can be at- 
tained without entanglement (with the EPR beams replaced 



by the vacuum). Braunstein et al. 16(] proposed a more rig- 
orous benchmark for CV teleportation. They considered the 
fidelity between an input state jV'in) and the corresponding 
state. Pout, outputted from a measure-and-prepare channel, 
i.e.. Pout = X]y P(xl'/'in)Px' whcrc p^ is the reconstructed 
or guessed state l3ased on outcome x of the measurement and 
p(xlV'in) is the conditional probability of obtaining x given 
that the signal state is iV'in)- Their benchmark is given by the 
maximum fidelity averaged over a conveniently chosen set fi 
of input states. It should be stressed that the choice of O plays 
a significant role. It should contain necessarily more than one 
state (otherwise the test is passed by a trivial classical strat- 
egy). The maximum fidelity tends to decrease with the size of 
the test family. When the input states are drawn from a whole 
d-dimensional Hilbert space, the optimal fidelity is known to 
be JF = 2/{d + 1) \M. In the case of CV systems, d ^^ oo, 
this would mean that any non-zero fidelity would certify that 
quantum resources are being used. However, it is utterly unre- 
alistic to assume that one can test the channel with such a large 
family of input states. In order to have realistic and practical 
thresholds, Braunstein et al. j^ chose Vt to be the set of coher- 
ent states with normally distributed amplitude around |q;| = 
and a fixed given variance, and gave a classical strategy that 
recently has been proved to be optimal by Hammerer etal. \^. 

Adesso and Chiribella [3|] have very recently proposed a 
quantum benchmark suited pure single-mode squeezed states 
and derived upper and lower bounds for squeezed thermal 
states. Their quantum benchmark is taken to be the maxi- 
mum averaged fidelity over all the outcomes x of an optimal 
measurement and over an ensemble Vt of input states pin {r) 
whose squeeziiig parameter r is distributed with a given prior 
probability (in 12[] a benchmark for pure squeezed states drawn 
from a microcanonical ensemble was considered). In this ap- 
proach, (i) the output state is assumed to belong to the input 
family, i.e., is of the form p{r^), and (ii) the "verifier" has to 
know what outcome x has been obtained at every single use 
of the hypothetical measure-and-prepare channel. This makes 
the so-defined fidelity difficult to compare with the fidelity of 
the real channel, which does not even have to involve mea- 
surements. As to (i), we will see that this choice of output 
state is typically sub-optimal and, hence, the coiTesponding 
fidelity might not be a strict quantum benchmark, as the au- 
thors noticed. 



Here we propose a benchmark based on a general phase- 
covariant family (pure or mixed): 
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(1) 



where i^(pi,P2) = (trVv^i^W^) = ^^AyfP\^fP2\) 
is the fidelity, p^ is the input test state defined as p^ = 
U{(t>)pQU{(i>)'^ with U{(t>) = E„e"^"|n)(n| (\n) are Fock 
states), and pav(0) = X]y ?'(x|P0)Px i^ '^^ ^'■^'■^ outputtedby 
the measure-and-prepare channel. In contrast to (i) above 121, 
[31], guessed states p^ are not constrained to belong to the 
input-state family. In contrast to (ii), pav is the output of the 
channel, with no reference to the way it has been implemented 
or to any further information. Since the fidelity for mixed 
states is non-linear Adesso and Chiribella's benchmark and 
ours will typically differ From the convexity of the fidelity 
it follows that theirs is a lowerbound to T and hence it is ar- 
guably not a proper quantum benchmark. A benchmark based 
on phase-covariant states is appealing from an experimental 
point of view, since phases are easy to vary without affecting 
other relevant parameters (e.g., in the presence of losses, vary- 
ing the degree of squeezing leads to a change in purity of the 
test states, which is very hard to compensate). 

Given a strategy characterized by POVM elements { O^ = 
l^xX^xl} (they can be taken to be of rank one without loss 
of generality and \^y^ are not normalized), and corresponding 
guesses {p^, one defines a covariant strategy by {O^^.b ~ 
l/(27r) UeO^Ul} and {px.e = Uep^Ul}. One can show that 
the optimal strategy can always be chosen to be covariant lll7ll 
and 



T={tl\^/p^y^\f; /9av=/ 
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Note that aside from the group parameter 6, one needs to 
specify also the "seed" for both POVM (O^) ^i^^ guess 
states ipx)- For single-seed strategies, the completeness re- 
lation fixes the POVM, which turns out to correspond to the 
phase-measurement ISl): |^) = ^^ |?i) (up to some arbitrary 
phases). 

The classical fidelity (|2]) can be conveniently written as 



T = max ( tr_B y trA VPo ® y/PoK^/lh) ® VPo 1 , (3) 

where tr^ (similarly tr^) stands for the obvious partial trace 
and K — Jdd ^ O^.e ® Px-0- Optimizing the classical strat- 
egy amounts to maximizing the trace squared in Q over the 
set of positive operators acting on TC ^Ti, that are separable, 
invariant under bilateral transformations Ug®Ug, and that ful- 
fill IybK = 1a- For pure states, Eq. ((Sj can be simplified to 
give J^ = (V'o|('0o|-f^|V'o)IV'o>, with po = l^oXV'ol- This 



leads to the pure-state estimation approach introduced in 1112 1. 
For a given POVM with seeds {|Cx)} the optimal fidelity can 
also be written as. 



Jc- = ^sup(?/'xl^xlV'x 
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where ^x = / d(/>/(27r)|(exlV'</,)PlV'</,)(V'^l- Thus, if the 
POVM is fixed, the maximum fidelity is given by the 
largest eigenvalue of the operator A^ and the optimal guess 
seed, |i/;^), by the corresponding normalized eigenvector llj]. 

If one restricts the guess-states to be in SI (as done in ISISti), 
things simplify considerably, specially for pure states where 
no optimization is required since the optimal POVM is known 
to be the phase-measurement 181]. In our case, no assumption 
about the optimal POVM nor about the guess is made and we 
have to resort to more powerful techniques. 

Semi-definite programming (SDP) |J9t| is an area of convex 
optimization that was developed in the last decade and that 
has recently found several applications in the field of Quan- 
tum Information ild.lll[ll2ll . Its aim is to minimize a linear 
objective function subject to semidefiniteness constraints in- 
volving symmetric matrices that are affine in the variables; 
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x,F, > 0, (5) 



where Ft are hermitian matrices of arbitrary dimension, and 
the inequality means that F{x) is positive semidefinite. There 
are a number of freely available software packages to solve 
SDP problems. In this work we have used the YALMIP mat- 
lab toolbox 1 13] together with the SDPT3 solver iH . 

For pure input states, the maximization in (O can be im- 
mediately cast into a SDP problem, by writing the matrices 
K and po in a basis of Hermitian matrices. Whereas the 
positivity and bilateral invariance constrains are already in 
the desired form, the separability condition need be imposed 
through a hierarchy of constrains based on PPT symmetric ex- 
tensions ilOtl . In this work we will stick to the lowest level of 
this hierarchy, imposing only positivity under partial transpo- 
sition (PPT), i.e., K^ > 0. Since positive partial transposi- 
tion provides a necessary, but in general not sufficient, condi- 
tion for separability, the resulting optimal value, T^ , gives an 
upper bound to T, and hence still provides a valid quantum 
benchmark. 

When the input test states are mixed, the objective func- 
tion in (O becomes truly non-linear and the optimization 
problem does not immediately fall into the SDP category. 
However, one can make use of Uhlmann's theorem 11511 



and recast J^ in ^ as J^ = max>p^^, |(*o|^ati)P = 
min„^^(- (4-0 1 CTav |*o)), where |*o) andffav = |*av) (*av| 
are purifications of po and pav respectively. Without loss of 
generality, the purity condition on cTav can be lifted. With 
this, the objective function becomes a linear function of the 
optimization variables. Constraints are also of the SDP form 
in ©: (i) trstTav == Pav = trA(po «) i K); (ii) CTav > 
and tr (Tav = 1; and (iii) the same conditions on K as above. 

Qubit states: It will become apparent as we proceed, that 
analytical solutions to the benchmark problem for general 
mixed states are exceedingly hard to obtain. A remarkable 
exception are qubit states, which we discuss next. 

The input-state family is defined by p^ = Uz,tj,poUl j,, 
where C/^ ^ is a rotation of an angle cj) around the z- 
axis (similarly for Uy,e, etc.), and the seed input state is 



Po = UyM{pn){]\ + {l-pmmulg, for fixed azimuthal 
angle and probability p = (l+77)/2 (so that r/ is the modulus 
of the Bloch vector). To calculate the quantum benchmark for 
this family of qubit states, we use ^. In the up/down basis the 
most general K can be written as i^ = blockdiag(a,i?,c), 
where i? is a 2 x 2 positive semi-definite matrix and a, c 
are non-negative numbers. For C^ (8) C^ the partial trans- 
position criteria provides a necessary and sufficient condi- 
tion for separability, i.e., K is separable <^ K^ > <=> 
ac — I-B12P > 0. Finally, the condition tr^A' = 1a im- 
plies a + Bii — c + B22 = 1. A tedious but straightforward 
calculation leads to the optimal K (that maximizes T): a = 
cos2((/2), c = siii^(C/2), B12 = B21 = y/ac, where the 
(azimuthal) angle ^ is 



7?2sin2 0+J{l- 772) (4 _ ,^2 sin2 61) 

C = arctan ^-- . (6) 

2t] cos 

It is a simple exercise to check that these values 
of a, B and c correspond to the phase-covariant POVM 
{l/TrUz,^Uy^:^\'\){'\\Uyj^Ul^}, with the associated guess 

10) = C/z^^C/y.^lt). The resulting benchmark for qubit States 
can be cast as 



T 



1 



77 



008 6* 

cos^ 



(7) 



Some remarks are in order; (i) for any outcome (f) the cor- 
responding guess is a pure state and, hence, (ii) it does not 
belong to the original family; moreover, (iii) its Bloch vector 
is not proportional to that of the signal state {( ^ 9, unless 
6 = 7r/2). In quantum state estimation it is usually assumed 
that the guess after a measurement is one of the possible in- 
put states. Our result show that this assumption is not always 
well-founded, (iv) Recalling the definition of K right after 
Eq. (O, we see that only one value of x (a single "seed") is 
required for qubits. Property (iv) is specific of C^. For C* 
one can already find examples where the phase-measurement 
is not optimal. 

We now move to the benchmarks for CV gaussian states. 
We consider displaced and squeezed thermal states that 
are obtained by the action of the displacement, D{a) = 
exp[a(a''' — a)], and the squeezing, S{r) = exp[|(a^ — a'^^)], 
operators over a thermal state pp = {l — c^^)e^^°' "• of purity 
/i = tanh/3/2 (with a > and r > 0). 

CV pure states: We start by computing the bound T^ us- 
ing the SDP approach for coherent states. For this purpose we 
use a truncated Fock basis and approximate low amplitude 

coherent states by |a) « e~" ^^ X]n=o '^"/'^'^l"')- ^^8" 
ure [T] (dots) shows the results for coherent states with mean 
photon number between zero and a^ = 10 (the truncation 
error within this range of a and A^ = 23 is always lower 
than 10^^). In addition we can calculate the optimal fi- 
delity when restricting to the single-seed covariant POVM. 
With this choice of POVM the problem reduces to calcu- 
lating the maximum eigenvalue of the matrix in ^, A ~ 



g-2|a| ^^«2„/^, i^^^^^^^j^ ^here |0„) = Y.t, V (?) 10- 
The values for different input intensities are shown in Fig. [T] 
(solid line) and agree with those obtained from the SDP opti- 
mization. This indicates that the benchmark given by T^ is at- 
tainable with the phase-measurement (at least within the pre- 
cision of our numerical analysis). We note that the eigenvector 
of A with largest eigenvalue, which is the optimal guess, re- 
sembles a coherent state but is strictly different. The dashed 
line in Fig. [T]is the fidelity obtained when the guess is forced 
to be a coherent state and, although it has a similar behavior, 
it shows a clear gap with the optimal bound. In this case no 
optimization is required and J- can be obtained by numerical 
integration of 



T 



/gi(eiac^^)n 



(8) 



From the above equation it is possible to find the analytical 
value of T in the limit a -^ 00. In this regime, the out- 
come probabilities can be approximated by | (CI "6**^)12 ~ 
■\/2a2/7rcxp[— 2a202j^ where we have used the Gaussian 
limiting expression of a Poisson distribution, Kalae**^)!^ « 
cxp[— a202j^ ^jjjj extend the range of integration to (—00, 00) 
to find Ta~,oo ~ ^2/3. 




FIG. 1: Plot of the average fidelity for pure coherent states (left) 
and pure squeezed states (right). Dots: bound obtained imposing 
K^ > (see main text). Sohd: phase-measurement and optimal 
guess. Dashed: phase-measurement and guess from input family. 
Dotted: extrapolation of dashed line to infinite squeezing, A = 1. 

The analytical value of JF in the large a limit can also be 
obtained in the general case of unrestricted guess states. In 
this case one has to calculate the maximum eigenvalue of A 
defined in Eq. Q, which can be done by calculating the limit 
p ^ 00 of its p-norm \\A\\p = (trAP)i/P. We have 

i\\A\\pr = tTAP= ff[dcj,,p{x\^,){a,\a,+^), (9) 

where a^+i = ai. Using the above approximation on the 
outcome probabilities and {ai\aj) « cxp{—a'^[i{(j>i — (j)j) + 
l/2{4>i - 0j)^]}, we obtain 



tiAP 



2a' 

TT 



2\P/2 



dP(j)C- 



-0 -Cp-tt) 



2P 



-y/det Cp 



(10) 



where Cp is a symmetric matrix with elements [Cp]i 



6S, 



^i+i.j 



Si,iSj^p for i < j. It is convenient to write 



Ct, 



Mr, 



2atp, with a = 3. Then, dct(Mp + 2at^ 



Qp{a) is a (characteristic) polynomial in a. It easy to check 
that Qp(a)_^2[Tp(a)- 1], where Tp(a) = [(a + ^/t?^^)P + 
[a — \/a? — 1)^] /2 are the Chebyshev polynomials of the first 
kind. In the limit p ^ oo the second term in Tp{a) becomes 
neghgible, hence Qp(3) « (3 + \/8)p = (\/2 - 1)^^, and 

T = lim 2(dct C„)^/2p ^ 2(\/2- 1) « 0.8284, (11) 

p — *oo 

which indicates that the difference between the fidelities cor- 
responding to restricted (guess in 51) and unrestricted (general 
guess) strategies persists also in the asymptotic regime. 

For pure squeezed states we proceed as be- 
fore by working in a truncated basis: |A) = 

(1 - ^')'/' ELo (-^/2)" ^fWY■|n\ \2n). Figure [U 
shows the SDP results, together with the bound obtained 
from the phase-measurement with optimal guess (the max- 
imum eigenvalue of A), and from the phase-measurement 
with (restricted) squeezed guess. The latter is obtained by 
numerically integrating 



T- 



d(j){l-X^) 



2^2 



1 + (i^^sin^.; 



E 



V(2^ /Ac^ 
nl \ 2 



(12) 



The SDP bound and the phase-measurement fidelity agree 
within numerical precision. As for coherent states, restricting 
the guess states be in il lowers the bound substantially. The 
latter bound can again be computed in the limit A ^ 1 by 
noticing that the dominant behavior of the sum in Eq. ( fT2] i is 
dictated by the large n terms 11611 . Using the Stirling's ap- 
proximation, the sum in ( fT2] i can be written as a polyloga- 
rithm function Lii/4(Ae^"^)/7r^'^. Taking into account that 
the dominant contribution to the integral comes from the re- 
gion ~ 0, we get JF ~ 0.58. 
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FIG. 2: Quantum benchmarks for displaced thermal states (left) and 
squeezed thermal states (right) for different purities: i-i — 1 (solid 
line), /i = 0.95 (dots), /i = 0.8 (dash) and ^ = 0.7 (dash-dot). 

CV mixed states: The case of mixed states is remarkably 
more complex and we have to entirely rely on numerical anal- 
ysis. Figure |2] shows the SDP results obtained for displaced 
thermal states and squeezed thermal states of different purity 
fi lllSfl . The truncation errors at the higher values of r and a^ 
are of the order of a few percent, but the displayed values still 
provide a good upper bound because truncation effects tend to 



lower the curves. We observe that decreasing the purity has 
the effect of increasing the fidelity. Thus our benchmark is 
specially suited for test states of moderate temperature. It is 
worth mentioning here that if the guess is restricted to belong 
to the input family (not shown in plot), the effect is the oppo- 
site: pure states provide higher fidelities than mixed states. 
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